Electronic structure and ionicity of actinide oxides from first principles calculations 
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The ground state electronic structures of the actinide oxides AO, A2O3 and AO2 (A=U, Np, 
Pu, Am, Cm, Bk, Cf) are determined from first-principles calculations, using the self-interaction 
corrected local spin-density (SIC-LSD) approximation. Emphasis is put on the degree of f-electron 
localization, which for AO2 and A2O3 is found to follow the stoichiometry, namely corresponding to 
A 4 + ions in the dioxide and A 3+ ions in the sesquioxides. In contrast, the A 2+ ionic configuration 
is not favorable in the monoxides, which therefore become metallic. The energetics of the oxidation 
and reduction of the actinide dioxides is discussed, and it is found that the dioxide is the most 
stable oxide for the actinides from Np onwards. Our study reveals a strong link between preferred 
oxidation number and degree of localization which is confirmed by comparing to the ground state 
configurations of the corresponding lanthanide oxides. The ionic nature of the actinide oxides 
emerges from the fact that only those compounds will form where the calculated ground state 
valency agrees with the nominal valency expected from a simple charge counting. 



I. INTRODUCTION 

Actinide oxides play a dominant role in the nuclear 
fuel cycle. [l[ For many years, uranium dioxide has been 
the main fuel component in commercial nuclear reactors. 
The " burning" of UO2 results in considerable amounts of 
Np and Pu isotopes, as well as smaller quantities of minor 
actinides such as Am, Cm, Bk and Cf. In original, "once 
through" reactors, the highly radioactive waste that is 
produced resulted in very troublesome long time storage 
requirements: particularly of Pu. However, it was soon 
realized that the Pu could be reprocessed from the spent 
fuel and used as alternative fuel in a new generation of re- 
actors. The environmental and energy production issues 
apart, the fact that Pu is obtained from the decommi- 
sioning of nuclear weapons is yet another important con- 
sideration concerning its use as nuclear fuel. A mixture of 
UO2 and Pu02, the so-called mixed oxide (MOX), where 
Pu is blended with either natural or depleted Uranium, 
constitutes the preferred, Pu containing, fuel in existing 
nuclear reactors. Lately, Np, as well as the minor ac- 
tinides that accumulate during nuclear reactor operation, 
are also being considered for reprocessing. [2j Once sep- 
arated from the spent fuel, these actinides can either be 
incorporated in durable ceramic waste for safe long time 
storage (and possible later recovery) or transmuted from 
long-lived isotopes to less radiotoxic short-lived isotopes 
through irradiation, thus taking part in the fuel cycle and 
reducing the long term nuclear waste management prob- 
lem. Again oxides are being considered both with respect 
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to the materials being used as fission/transmutation tar- 
gets (A0 2 ) 1 and for the direct storage in the shape 
of durable ceramic glasses (A0 2 , A 2 3 ). For exam- 
ple PuC>2 has long been the compound of choice for de- 
positing Pu in long time repositories, given the observed 
stability with respect to oxidation. [J] 

Regarding their behaviour under both reactor opera- 
tion and storage conditions, it is crucial to understand 
the thermochemistry, thermophysics, and materials sci- 
ence of the actinide compounds. Given the toxicity of 
the materials involved, computer simulations, such as 
thermodynamic modelling [5[ or, as in the present pa- 
per, electronic structure calculations can provide funda- 
mental insights at a level not achieveable through exper- 
iment alone. Here we wish to focus specifically on the 
/-electrons, their contribution to the groundstate elec- 
tronic properties of the actinide oxides, the role they play 
with regards to stability towards oxidation, and their be- 
haviour under ionic bonding conditions. 

When modelling the electronic structure of actinide 
materials, the most distinguishing feature is the increas- 
ing importance of correlations across the series from U 
to Cf, as the nature of the /-electrons changes from de- 
localized in the early actinides to localized in the later 
actinides. @, 0] Electronic structure calculations, based 
on the local spin density approximation (LSD A), do not 
take into account strong on-site correlations beyond the 
homogeneous electron gas, and therefore can not ade- 
quately describe the localized phase of actinide materi- 
als. Thus the LSDA, Q or even the generalized gradi- 
ent approximation (GGA) Q (which extends beyond the 
LSDA by taking into account charge density gradients), 
wrongly predicts a metallic ground state for UO2, Pu02 
and PU2O3, although the equilibrium lattice parameters 
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and cohesive properties are found to be in rather good 
agreement with experiment. 0, [l(| [H| 

A number of schemes have been developed that aug- 
ment the standard bandstructurc framework to include 
the effects of strong correlations on the electronic struc- 
ture. In the LDA+U approach [l^] an effective Coulomb 
parameter U is introduced that separates the /-manifold 
into the upper and lower Hubbard bands and removes 
/-degrees of freedom from the Fermi level. The more ad- 
vanced dynamical mean-field theory (DMFT) approach 
provides a description of the competing trends towards 
localization on the one hand and itineracy on the other 
hand by taking into account the local quantum fluctua- 
tions missing in the static LDA+U treatment [HI, [l4| ; 
albeit still at the cost of the introduction of the U- 
parameter. The hybrid density functional [l5| theory 
implements an exchange-correlation functional where a 
fraction of the exact non-local exchange interaction from 
Hartree-Fock theory is mixed with the local or semi-local 
exchange energy of LSDA or GGA with the result that 
the troublesome effects of the known self-interaction er- 
ror present in the standard LSDA and GGA calculations 
are reduced. 

The self-interaction corrected (SIC)-LSD approach [l6[ 
used in the current work removes the self-interaction er- 
ror that occurs in the LSDA, thereby leading to an im- 
proved description of the static Coulomb correlations of 
the f-electrons. The self-interaction correction associates 
an energy gain with electron localization, which com- 
petes with the opposing trend of band formation, pro- 
viding a dual picture of combined localized and band 
like /-electrons. The method is fully ab-initio as both 
kinds of electrons are treated on an equal footing, with 
no adjustable parameters. A comparative study of MnO, 
involving SIC-LSD, LDA+U, and the hybrid functional 
methods, was published by Kasinathan et at. [I?]] 

The present paper is organized as follows. In the fol- 
lowing section, we give a short introduction to the SIC- 
LSD band structure method. In section III, we present 
our SIC-LSD results for the ground state properties of 
(A) the monoxides, (B) the sesquioxides, and (C) the 
dioxides where we also consider oxidation/reduction en- 
ergies. In section IV, we give a summarizing discussion 
of the results, also concerning the relation between /- 
electron localization and oxidation by comparing to the 
lanthanide oxides. The conclusion of our paper is pre- 
sented in section V. 



II. THE SIC-LSD METHODOLOGY 

The LSD approximation to the exchange and corre- 
lation energy introduces an unphysical interaction of an 
electron with itself [l6[ which, although insignificant for 
extended band states, may lead to uncontrollable er- 
rors in the description of atomic-like localized states, 
for example the /-electrons in the later actinides. The 
SIC-LSD method [H, Oil corrects for this spurious self- 



interaction by adding to the LSD total energy functional 
an explicit energy contribution for an electron to localize. 
The resulting, orbital dependent, SIC-LSD total energy 
functional has the form 
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Here the sums run over all occupied electron states ip a . 
As usual, the LSD total energy functional (2) is decom- 
posed into the kinetic energy, the Hartree energy, the in- 
teraction energy with the atomic ions, and the exchange 
and correlation energy. The spin-orbit interaction (3) 
couples the band Hamiltonian for the spin-up and spin- 
down channels, i.e. a double secular problem must be 
solved. The spin-orbit parameter, 
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in atomic Rydberg units, is calculated from the self- 
consistent potential. The self-interaction energy (4) con- 
sists of the self-Coulomb and self-exchange-correlation 
energies of the occupied orbitals ip a with the orbital 
charge density n a . 

For itinerant states, the self- interaction S^ IC vanishes 
identically, while for localized (atomic-like) states 8^ IC 
may be appreciable. Thus, the self-interaction correc- 
tion constitutes a negative energy contribution gained by 
an electron upon localization, which competes with the 
band formation energy gained by the electron if allowed 
to delocalize and hybridize with the available conduction 
states. Different localized/delocalized configurations are 
realized by assuming different numbers of localized states 
- here /-states on actinide-atom sites. For s- and p-states, 
(5f IC is never competitive compared to the correspond- 
ing gain in band formation energy, and turns out to be 
positive. Since the different localization scenarios con- 
stitute distinct local minima of the same energy func- 
tional, E SIC ~ LSD , their total energies may be compared 
and the global energy minimum then defines the ground 
state total energy and the valence configuration of the 
actinide-ion. This latter is defined as the integer number 
of electrons available for band formation, 
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atomic core electrons, and N$ic is the number of SIC- 
localized /-electrons. In the remainder of this paper we 
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will be using two interchangeable nomenclatures, f n and 
A m+ , to describe the configuration of the actinide ion, 
implying n = Nsic an d m = N va i, respectively. The to- 
tal number of /-electrons may be larger than n, since, in 
addition to the n localized /-states, the band states also 
contribute to the total /-count on a given ion. Note that 
our calculated valencies refer to the number of actinide 
electrons that contribute to bonding, and thus do not 
necessarily coincide with the nominal (ionic) valency of a 
compound: For PuC>2 for example, the Pu 4+ would agree 
with an ionic picture, while the Pu 3+ , and Pu 5+ valency 
configurations would indicate more covalent behaviour. 

The SIC-LSD approach is fully ab-initio, as both lo- 
calized and delocalized states are expanded in the same 
set of basis functions and are, thus, treated on an equal 
footing. If no localized states are assumed, e sic ~ lsd 
coincides with the conventional LSD functional, i.e., the 
Kohn-Sham minimum of the E LSD functional is also a 
local minimum of E SIC ~ LSD . 

Given the total energy functional ]t;SIC~lsd ^ ^ c 
computational procedure is as for the LSD case, i.e. 
minimization is accomplished by iteration until self- 
consistency. In the present work, the electron wave- 
functions are expanded in the linear-muffin-tin-orbital 
(LMTO) basis functions. [20| The atomic spheres approx- 
imation (ASA) is used, whereby the crystal volume is 
divided into slightly overlapping atom-centered spheres 
of a total volume equal to the actual volume. A known 
shortcoming of the ASA is that different crystal struc- 
tures have different degrees of overlap of the ASA spheres 
resulting in substantial relative errors in the evaluation 
of the total energy. While this inhibits the comparison 
of energies of different crystal structures, when compar- 
ing the energies of different localization scenarios within 
the same crystal structure, the ASA error is of minor in- 
fluence. To improve the packing of the structure empty 
spheres have been introduced on high symmetry intersti- 
tial sites. 



III. RESULTS 



A. Actinide Monoxides 



Background information 



There exists to date no convincing evidence that ac- 
tinide oxides can form in the 1:1 stoichiometry. The ex- 
perimental lattice parameters that we cite in Tablc[T]comc 
from early reports on these compounds, and have so far 
not been confirmed. There have been no claims of bulk 
UO having ever been synthesized, and reports of an UO 
surface phase on U metal for low exposures to O . l2ll. I22I 
and UO thin films [23| could not be reproduced. 24 1 It 
has been suggested that the observed thin films actually 
represent uranium oxynitrides (UN^Oi-^) and oxycar- 
bides (UCzOi-z), that form in the presence of N2 or C, 
and at low oxygen pressure. [25| It has similarly been con- 



TABLE I: Actinide monoxide data. Column 2: Groundstate 
configuration. Column 3: Energy difference in eV between 
the ground state and the ideal ionic divalent configuration. 
Column 4: Groundstate density of states at the Fermi level 
(in units of states per eV and formula unit). Column 5: Cal- 
culated lattice parameters aQ alc (in A) . Column 6: Experi- 
mental lattice parameter ag Xp (in A), where known (measure- 
ments on UO, and NpO refer to thin film data, measurements 
on PuO, AmO. and BkO refer to bulk data). 



AO 


Config. 


Egs-E// 


n(E F ) 


calc 

a 


exp 

a 


UO 


f 1 (U b+ ) 


-1.93 


5.8 


4.94 


4.92° 


NpO 


f a (Np 4+ ) 


-1.73 


1.1 


4.99 


5.01" 


PuO 


f b (Pu 3+ ) 


-0.58 


5.9 


5.13 


4.960 c 


AmO 


f° (Am a+ ) 


-0.14 


7.3 


5.14 


5.045 c 


CmO 


{' (Cm J +) 


-1.14 


2.6 


5.01 




BkO 


f (Bk a +) 


-0.65 


1.4 


4.97 


4.964 d 


CfO 


f (Cf + ) 


-0.20 


9.5 


4.97 




EsO 


f iu (Es J+ ) 


0.00 


4.0 


4.92 




EsO 


f u (Es 2 +) 


0.00 


0.0 


5.06 





'Reference 23 b Rcfcrcncc 30 Reference 28 Reference 29 



eluded that neither bulk NpO, nor a corresponding NpO 
surface phase will form. [20, [27[ Pre par ation of PuO and 
AmO, [28[ as well as possibly BkO, |29| has been neither 
substantiated, nor has it been dismissed. 



2. SIC-LSD electronic structure 

We have calculated the electronic structure of the 
monoxides with the SIC-LSD method in order to estab- 
lish the ground state properties for the hypothetical NaCl 
structure. Ferromagnetic arrangement of the spins has 
been assumed in these calculations. The results are sum- 
marized in Table HI We find the trivalent configuration 
to be energetically most favourable for all the monoxides, 
except UO and NpO that respectively prefer the U 5+ 
and Np 4+ ground state configurations. Concerning UO, 
it should be noted here that even though we find the U 
f 1 configuration to be the ground state, this energy sce- 
nario is close to degenerate with the fully itinerant LSD 
(/°) and f 2 scenarios, as shown in figure [TJ This indi- 
cates that the /-electron manifold lingers between the 
localized and delocalized pictures, i.e. correlations are 
strong, but not to the degree for full localization to occur. 
This agrees with earlier results by Brooks et aZ.,[3l[ where 
the electronic structure of UO was calculated assuming 
itinerant /-electrons, and which resulted in a calculated 
lattice parameter of 4.88 A, i.e. only slighty overbinding 
with respect to the "experimental" value of 4.92 A. The 
SIC-LSD calculated lattice parameters of the f 1 and / 2 
configurations, respectively 4.94 A and 4.99 A, indicate 
a slight overlocalization. Overall, the calculated lattice 
parameters are in rather good agreement with the early 
experimental values. However, as there is no convincing 
subsequent experimental evidence in support of the fact 
that monoxides really exist in nature, one should not put 
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4.8 5 5.2 

Lattice parameter (A) 



FIG. 1: Total energy of UO in the rocksalt structure assuming 
different localized/delocalized /-electron scenarios. 




Energy (eV) 



FIG. 2: Density of states for CfO in the ground state config- 
uration with nine /-electrons localized on the Cf 3+ ion. The 
majority (minority) spin components are displayed as positive 
(negative) values. The energy is measured in eV relative to 
the Fermi level, which falls amidst the minority Cf / bands. 
The localized Cf / states are not shown. 



too much weight on the agreement between theory and 
experiment. 

For the monoxides beyond UO, the LSD configuration 
never becomes even remotely energetically favourable. 
For NpO a tetravalent ground state (Np 4+ = Np(/ 3 )) is 
found, whereas a trivalcnt ground state configuration is 
established for the remaining monoxides. In their respec- 
tive ground states, the monoxides are metallic, as can be 



seen from the non-zero density of states (DOS) at the 
Fermi energy (column 4 of Table QJ. As a representa- 
tive example, the DOS of CfO is depicted in Fig. [5] 
The O atom has two unoccupied p-states, and in the cor- 
responding monoxide, the p-band can accomodate two 
electrons from the actinide atom through charge trans- 
fer and hybridization, whilst the remaining valence elec- 
trons, including the delocalized /-electrons, start filling 
the conduction band, with the Fermi level pinned to the 
narrow /-band. The ionic insulating picture would be 
realized if additional /-states preferred to localize, i.e. in 
the divalent configuration. In Table [T] (column 3) the cal- 
culated energy differences, Egs-E//, between the ground 
state configuration and the divalent configuration arc 
shown. It is clear that the nominal ionic A 2+ 2_ scenario 
does not become energetically favourable for any of the 
monoxides. Nevertheless this energy difference is seen 
to decrease from UO to AmO, and again from CmO to 
CfO, which is in agreement with the trends towards in- 
creasing localization due to the actinide contraction. At 
AmO the trivalent configuration is thus only marginally 
more favourable than the divalent half-filled shell con- 
figuration, and at EsO degeneracy occurs between the 
two configurations. The overall non-ionic ground state 
picture for the monoxides that emerges from our SIC- 
LSD calculations confirms the results from early molecu- 
lar cluster calculations on the heavy actinide monoxides 
by Gubanov et aL,[32j where considerable covalency due 
to mixing of O-p and A-/ orbitals was found. 



B. Actinide Sesquioxides 

1. Background information 

Bulk phases of U2O3 as well as NP2O3 do not exist in 
nature, and have never been synthesized, but thin films 
of NP2O3 have been found to form on the surface of Np 
metal. [271] From Pu onwards, the sesquioxides are stable, 
and have been synthesized. Beyond Pu, the sesquiox- 
ides crystallize in three different crystal structures, re- 
spectively the hexagonal La203 structure (A-form), the 
cubic Mn 2 03 structure (C-form), and the monoclinic 
Sm 2 03 structure (B-form). Pu 2 03 has been synthe- 
sized only in the A- and C-forms. The XPS measure- 
ments on sesquioxides from PU2O3 to Cf 2 03 have been 
reported, [H, 0, [H] and the absence of features at the 
Fermi level points towards the localized nature of the 5/ 
electrons in these compounds. This indicates that they 
are semiconductors or insulators, in agreement with the 
ideal nominal picture of A^O 2- , although no values for 
the energy gaps can be found in the literature. 

There exist relatively few calculations of the elec- 
tronic structure of the actinide sesquioxides. Prodan 
et aZ.[36[ have studied PU2O3 (A-form) with the help 
of hybrid density functional theory, comparing a num- 
ber of functionals. They have found that unlike in the 
LSD and GGA approximations, using the Hcyd-Scuseria- 



TABLE II: Actinide sesquioxide data: Column 2: Ground 
state configuration of actinide-ion. Column 3: Calculated 
energy gap, E gap , (in e V) . Columns 4 and 5: Calculated, 
Vcaic and experimental [30(1 . V e x P , equilibrium volume(s) in 
in units of A 3 per formula unit. 



Compound | Ground state | E gap | V ca ; c | "Vexp 


C-type sesquioxides 


U2O3 


r (u 4+ ) 


0.00 


83.17 




Np 2 3 


f J (Np 4+ ) 


0.00 


84.40 




Pu 2 3 


f° (Pu d+ ) 


0.75 


89.42 


82.73 


Am 2 3 


f (Am 3+ ) 


0.44 


88.54 


83.64 


Cm 2 3 


f (Cm J+ ) 


0.32 


86.98 


83.10 


Bk 2 3 


f (Bk d+ ) 


0.38 


83.41 


80.63 


Cf 2 Q 3 


P (Cf 3+ j 


0.47 


82.60 


79.59 


A-type sesquioxides 


Pu 2 3 


f° (Pu d+ ) 


2.43 


74.06 


75.49 


Am 2 3 


f (Am 3+ ) 


2.54 


73.34 


74.73 


Cm 2 3 


f (Cm J+ ) 


3.07 


72.40 


74.53 


Bk 2 3 


f 8 (Bk 3+ ) 


2.73 


70.10 


72.71 


Cf 2 Q 3 


f (Cf 3+ ) 


1.78 


69.33 


71.43 



Ernzerhof (HSE) screened Coulomb functional leads to 
an insulating anti-ferromagnetic solution, in good agree- 
ment with experiment, with the calculated gap value of 
2.78 eV (3.50 eV for the Perdew-Burke-Ernzerhof (PBE0) 
hybrid functional). The antiferromagnetic insulating na- 
ture of PU2O3 has similarly been retrieved from electronic 
structure calculations by respectively Jomard et aLpj| 
and Sun et al.^^, based on the (LDA/GGA)+U approx- 
imation. Their gap values are found to depend strongly 
on the value chosen for the electron-electron interaction 
U. For U=4 eV, and depending on the details of the func- 
tional, energy gaps ranging from 1 to 2 eV were obtained. 



2. SIC-LSD electronic structure 

In the present work the electronic structures of the 
sesquioxides from U2O3 to Cf2C>3 has been calculated for 
both the cubic C-form and the hexagonal A-form. The 
results are summarized in Table [XXJ The hexagonal A- 
type structure has space group P3ml (no. 164). The 
unit cell contains one formula unit with one Oxygen at 
the origin, two Oxygens at ±(1/3, 2/3, zo) and two ac- 
tinide atoms at ±(1/3, 2/3, 2, a)- While the c/a ratio has 
been measured for all the actindes from Pu to Cf, the 
internal parameters are only known for Pu (zo = 0.6451 
and zp u = 0.2 10 s * ;. 3(1 Hence we have performed our 
calculations for this structure, with the experimental c/a 
ratios but using the internal parameters of PU2O3 for all 
the other actinide sesquioxides. [38| The C-type structure, 
also known as bixbyite, has been approximated by the flu- 
orite AO2 structure, with 1/4 of the O atoms removed, 
i.e. from four formula units in a conventional simple cu- 
bic supercell the oxygen atoms at the origin and the cube 
center were replaced by empty spheres. No relaxation of 
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E-E . (eV) 
fermi v ' 

FIG. 3: Density of states for A-type Pu 2 3 in the trivalent 
ground state configuration. The energy is in eV relative to 
the Fermi level, which is situated midgap. The DOS is in 
units of eV -1 per formula unit. Upper panel: Total (black), 
Pu (violet), and O (light blue) DOS's. Lower panel: The 
corresponding Pu-/ majority (red), and minority (blue) spin 
decomposed DOS's. Only the itinerant /-states are shown. 



the atomic coordinates has been attempted. The triva- 
lent actinide configuration is found to be energetically 
favourable for all the compounds, apart from U2O3 and 
NP2O3, for which the tetravalent ground state configu- 
ration is preferred. Incidentally, these are also the only 
sesquioxides that do not occur in nature. In the trivalent 
ground state, the sesquioxides are found to be insulators, 
with energy gaps of around 0.5 eV for the C-form, and 
around 2.7 eV for the A-form. The corresponding DOS 
for the A-type PU2O3 is shown in Fig. [3l 

Only the ferromagnetic ordering was investigated for 
the cubic structure, whilst only the antiferromagnetic or- 
dering of spins was studied for the hexagonal structures, 
and the values for the insulating gaps in Table [Hi refer to 
these respective magnetic orderings. The bulk moduli for 
A-type Pu 2 3 , Am 2 3 , Cm 2 3 , Bk 2 3 , and Cf 2 3 are 
similar in magnitude, calculated to be respectively, 158, 
158, 168, 166, and 174 GPa. No experimental measure- 
ments of the bulk moduli seem to exist, but the value for 
PU2O3 lies within the range of values obtained by Pro- 
dan et al. [36[ and Jomard et al. [i"lT | for a suite of different 
functional (from 110 GPa to 181 GPa). As can be seen 
from Table [TTJ the calculated equilibrium volumes of the 
A-type sesquioxides are found to be in good agreement 
with experiment. This agreement is somewhat less sat- 
isfying for the C-type sesquioxides, where the volume is 
overstimated by up to 4 %, which might be related to our 
approximate treatment of the actual bixbyite structure. 
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C. Actinide Dioxides 

1. Background information 

The actinide dioxides from U to Cf have all been syn- 
thesized, crystallizing in the fluorite (CaF 2 ) structure. 
They are also the most relevant systems for applications 
and, as a consequence, are the most studied actinide ox- 
ide compounds, both exp er iment ally and theoretically. 

Jji EH EE El EE El, E3, EE S3, S3, El ES EE EE Si, El, 

l50j On the experimental side, a comprehensive summary 
of their electronic, magnetic, transport and optical prop- 
erties can be found in Troc et aLj48[ The absence of fea- 
tures at the Fermi level in the observed XPS spectra[35j 
indicates that all the dioxides are semiconductors or in- 
sulators. However, specific information regarding the gap 
is only known for UO2, NpC>2, and PuC>2- A combined 
XPS and BIS spectrum for U0 2 finds an 0(2p) -> U(6d) 
gap of about 5 eV.[49[ This value for the energy gap is 
similar to the one observed for TI1O2, but with the dif- 
ference that in UO2 two occupied rather well localized 
/-states are situated in the gap, with a 5/ 2 — > 5/ 1 6sd 
transition energy of 2-3 eV.|43,l51| As the nuclear charge 
increases from Th to Pu, the /-states move to lower en- 
ergies. In TI1O2 the empty /-states are situated in the 
Th-s<i derived conduction band. In UO2 and NpC>2 the 
/-states arc occupied and situated in the energy band 
gap, whereas in PuC>2 they are situated at the top of the 
O-p derived valence band.[33[ The measured activation 
energies, E Q =0.2 eV in U0 2 , E a =0.4 eV in Np0 2 , and 
E a = 1.8 cV in P11O2, confirm this trend. [50j 

Electronic structure calculations, in particular for UO2 
and P11O2, have shown that the itinerant /-electron pic- 
ture is not adquate for describing these compounds. It 
has emerged that the LSD approximation wrongly pre- 
dicts metallic behaviour, [1 [ll], El] demonstrating the 
need to go beyond the homogeneous electron gas in de- 
scribing the strong on-site /-electron correlations. So 
far, a number of calculations have shown that an insu- 
lating solution can be obtained when using improved de- 
scriptions of electron correlations such as those contained 
in the LDA(GGA)+U[ll|, EE ES El, DMFTE1, SIC- 
LSD[43[ and hybrid functional methodologies. [44j 

2. SIC-LSD electronic structure 

In the present paper the electronic structures of the 
dioxides from UO2 to CTO2 has been calculated using 
the SIC-LSD method. As can be seen from Table Mil a 
tetravalent ground state configuration is found for all the 
dioxides, except UO2 where the pentavalent (U(/ 1 )) con- 
figuration is energetically more favourable. Actually, the 
pentavalent U(/ 1 ) and tetravalent U(/ 2 ) configurations 
are energetically close, with the f 1 configuration being 
more favourable by some 100 meV. Our calculations re- 
fer to T=0 K, whereas the experimental evidence, which 
clearly indicates a tetravalent (insulating) UO2, mostly 



refers to room temperature conditions. Experiments also 
indicate a lattice expansion with temperature . l53l| but 
the f 1 to / 2 localization transition that our calculations 
seem to predict has not been observed experimentally. 
The prediction of an U-/ 1 goundstate configuration could 
possibly be related to the tetrad effect (multiplet forma- 
tion energy), EE EH which is ignored in our calcu- 
lation, and which favours the / 2 configuration over the 
f 1 configuration. In the following, when comparing to 
experiment, we will be referring to the tetravalent UO2 
configuration as the ground state configuration. 

All the dioxides (including the tetravalent UO2) are 
predicted to be insulators. The AFII magnetic order- 
ing (i.e. fcrromagnctically ordered planes stacked anti- 
ferrormagnetically along the [111] crystal direction) has 
been assumed in the dioxide calculations to which the 
SIC-LSD data in Table HH refer. With respect to the 
band gap of UO2, a value of E gop =2.6 eV is found, which 
is smaller than the experimentally observed p— >d gap, 
and which should not be compared to the optical gap 
referred to in the GGA+U (1.8 cV), and HSE (2.4 eV) 
calculations. The SIC-LSD, being a one-electron ground 
state theory, does not give accurate removal energies of 
localized states, on account of the missing screening and 
relaxation effects, predicting them at too high binding 
energies, in disagreement with spectroscopy. As a re- 
sult, the SIC-LSD calulations do not reproduce the cor- 
rect position of the occupied /-peak which is expected 
to be situated in the gap for both UO2 or Np02- The 
Mott-insulating character (/— >/ transition) of the early 
dioxides [40l l44j is thus not reproduced in our calcula- 
tions, which instead describe the entire series as charge 
transfer insulators (p—>f transition). From Pu02 on- 
wards, the occupied /-states are situated at or below the 
valence band maximum, and the charge transfer picture 
becomes adequate for describing the nature of the gap. 
There exist only a few calculations dealing with the elec- 
tronic structure of the heavy actinide dioxides. [HI, El, SE 
In the following we make a detailed comparison with the 
recent results obtained by Prodan et al. [44| using the hy- 
brid functional theory. 

Both the SIC-LSD and the hybrid functional calcula- 
tions determine the antifcrromagnctic ground state con- 
figuration as energetically most favourable. As can be 
seen in Table Hill our calculated energy gaps are smaller 
than the corresponding HSE gaps, but the trends agree, 
i.e. the gaps decrease from Pu02 to Cm02, and then 
again from Bk02 onwards. This reflects the gradual pro- 
gression towards lower energy of the unoccupied /-states, 
with increasing atomic number. The calculated lattice 
parameters arc in very good agreement with experiment, 
within 1 %, as seen in both Table HTT1 and figured) We 
should note here that the experimentally observed lattice 
parameters refer to room temperature measurements . f30j 
and are thus on average larger than the values one would 
expect for T=0 K by approximately 0.011 A. (59| The con- 
sistent ovcrcstimation of the lattice parameters is quite 
common for SIC-LSD, which has tendency to slightly 
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TABLE III: Actinide dioxide data: Column 2: ground state 
configuration. Columns 3 and 4: Energy gaps, E^F (this 
work) and E^,p B (reference 0) , respectively (in units of eV). 
Columns 5 and 6: calculated, slq IC , and experimental a,Q Xp 
lattice parameters (in A). Column 7 and 8: calculated, Bq c \ 
and experimental, Bo IC , bulk modulii (in units of GPa). 



the ground state. [60[ It has also been conjectured that 
a more complex ground state which includes some ad- 
mixture of the excited J = 1 state of Cm 4+ might be 
required. [64 
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ovcrlocalizc. 

The overall experimental situation with regard to 
CmO 2 is still not fully understood. [60j It has been 
shown that, for a range of off-stoichiomctric compounds 
Cm02_i, both the lattice parameter and magnetic mo- 
ment increase with x.[6ll.[63.[63j Given the expected mo- 
ments of Cm 3 + (5/ 7 ; /x e// =7.94 (j, b ) and Cm 4 + (5/ 6 ; 
Me//— /is) it was noted that this trend indicates an 
increase in the number of Cm 3+ impurities as we move 
away from stoichiomctry. Based on the ionic picture, one 
accordingly expects the Cm 4+ ground state configuration 
for the stoichiometric compound, which however seems to 
be contradicted by the fact that in samples very close to 
stoichimetry, [64| susceptibility measurements give the ef- 
fective moment as high as fi e ff= 3.36 [Ib- Attempts to 
explain this discrepancy, range from a possible existence 
of an impurity phase, not registered in the experiments, 
to assuming a covalent, rather than an ionic, picture for 



The lattice parameter measurements yield two differ- 
ent values for Cm02, depending on whether the short- 
lived curium isotope of mass number 244 is used (a = 
5.372 A),[65[ or the long-lived curium isotope of mass 
number 248 is used (a = 5.359 A).(H2] The reason for 
this discrepancy has been traced to the self-radiation in- 
duced expansion that occurs with the considerably more 
active 244 Cm isotope. [H H3| It is the 248 Cm0 2 lattice 
parameter, i.e. without the effect of self-heating, that 
should be of relevance when comparing to the calculated 
lattice parameters. 

From the SIC-LSD calculations we find the ground 
state of C111O2 to be tetravalent, and the corresponding 
lattice parameter is in good agreement with experiment. 
In figure [5] we show the DOS of Cm0 2 in the Cm 4+ (/ 6 ) 
ground state configuration. The compound is found to be 
insulating, with the Fermi level positioned between the 
completely filled O-p band, and the one remaining de- 
localizcd empty majority /-spin state that strongly hy- 
bridizes with the O-p states. The resulting electronic 
structure is thus quite different from the hybrid func- 
tional picture, where the HSE applied by Prodan et al. 
results in metallic Cm02- The authors suggest a cova- 
lent picture with significant Cm 3+ character, as a result 
of the Cm trying to achieve the stable half-filled /-shell 
configuration. Consequently some of the O p-states are 
charge transferred to the Cm-ion, resulting in the Fermi 
level cutting accross the top of the hybridized Op- Cm 
/ band. However, unlike for all the other actinide diox- 
ides, the lattice parameter of C111O2, calculated using the 
HSE functional, deviates considerably from the measured 
value, as can be seen from Fig. 2] 

With respect to the SIC-LSD calculations, one could 
envisage that a Cm 3+ configuration in Cm0 2 could be 
energetically favourable, given the associated half-filled 
/-shell. Due to one less /-state taking part in bonding, 
as compared to the Cm 4+ , this in return could explain 
the slight increase in the experimentally observed lat- 
tice parameter of Cm02- As it turns out, however, the 
total energy of the trivalent configuration is higher by 
0.68 eV, indicating that the Cm 4+ configuration is en- 
ergetically very stable. The reason for this stability is 
related to the fact that localizing an additional electron 
in Cm0 2 results in the Fermi level moving down into 
the p-band. This implies a depopulation of the p-band, 
as charge is transferred to the low lying /-levels, and 
the associated loss in respectively the Madelung and hy- 
bridization energies is significantly larger than the gain 
in the /-localization energy. This charge transfer picture 
is similar to the picture found in the hybrid functional 
theory. However, in variance to the conclusion of that 
work, we find the corresponding trivalent configuration 
to be energetically unfavourable. 

Although our calculated ground state of C111O2 finds a 
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localized / 6 Cm ion, the system turns out to be mag- 
netic The ideal / 6 ion has J = and hence also 

< S z >= and < L z >= 0, but this state is a linear 
combination of several 6-electron Slater determinants, 
which we cannot treat in our SIC-LSD scheme. Rather, it 
may be represented by a single Slater determinant hav- 
ing < S z >= 3 and < L z >= —3 (simulating antipar- 
allel S = 3 and L = 3), which constitutes the start- 
ing point for our calculations before iteration to self- 
consistency. Spin-orbit coupling and hybridization sig- 
nificantly distort this initial configuration, in fact almost 
quenching the orbital moment. The total spin and or- 
bital moment projections along the z-axis are calculated 
to be < S z >= 2.73 and < L z >= -0.25 (in units 
of h). This comes from contributions from the local- 
ized / 6 ion ( < S z >= 2.76 and < L z > = -1.15) and 
from the delocalized f-band states ( < S z >= —0.03 and 

< L z >= +0.93). The latter are largely tails of the 
O p-band states, which inside the atomic sphere around 
Cm attain /-character. The total number of delocalized 
/-character states is 0.86, leading to a total /-count of 
6.86, which reflects a significant /-hybridization with the 
O p bands. Assuming now, that the magnetic moment 
is given as fi =< 2S Z + L z > /is, we arrive at a to- 
tal moment of [i = 5.21/x_b per Cm atom, which exceeds 
the experimental moment of [J, e ff = 3.36/xs,:64] how- 
ever demonstrating that the insulating state of C111O2 
may be magnetic. The quoted experimental moment is 
extracted from the large temperature magnetic suscep- 
tibility, while the present theory is valid only at T = 0. 
Furthermore, as stated, the fact that quantum fluctua- 
tions (more than one Slater determinant representing the 
localized / 6 shell), are not possible within the SIC-LSD 
approach, might be a serious limitation with respect to a 



proper description of magnetic properties. Thermal fluc- 
tuations can be considered using the so-called local SIC 
(LSIC) approach, implemented in the multiple scattering 
theory, in combination with the coherent potential ap- 
proximation (CPA) and disordered local moment (DLM) 
theory. [H [63] 



3. Oxidation and Reduction Energies 

Uranium metal readily oxidizes to form a range of 
super-stoichiomctric oxide compounds, \JO-2+ x , as well as 
stoichiometric compounds, UO2, U3O8, UO3, and U2O5. 
Sub-stoichiomctric UO2-1 exists, but U2O3 for example 
is not found in the literature. In this section we will inves- 
tigate oxidation/reduction of the actinide dioxides, based 
on the SIC-LSD total energies involved in the different 
delocalization/localization transitions of the /-electrons. 
To model the oxidation process 



A0 2 + ±O n 



AO 



2.25, 



(6) 



one additional O is introduced into a four formula units 
AO2 supercell, i.e. we define the oxidation energy 



E ox = _ 



E(A 4 9 ) - £(A 4 8 ) - i-B(0 2 ) 



(7) 



The CaF2 structure is assummed to remain undistorted, 
i.e. relaxation effects are not taken into account. For the 
AO2 reduction process we consider the reaction 



AO, 



AOr. £ 



(8) 



In the supercell used for modelling the reaction, two O 
atoms are removed from the four formula units AO2 su- 
percell. Correspondingly, the reduction energy is 

E red = 1 [£(A 4 6 ) + E(0 2 ) - E(A 4 8 )} . (9) 

The compound A4O6 is basically a two formula unit 
supercell of the sesquioxide A2O3 in the previously de- 
scribed approximation to the bixbyite structure. 

The supercell total energies E(A40s) and E(A40g) 
have been evaluated at the ground state valency con- 
figurations established in the preceding sections for 
respectively AO2 and A2O3. For A4O9, the SIC- 
LSD calculations find the following ground state con- 
figurations: U^ + 9 , Np^ + 9 , Pu^ + Pu^ + O g , Am^+Og, 
Cmg+Cn^+Og/Cm^+Og, Bkjj+Og, and Cf^+Og. The 
corresponding total energy minima E(A40g) are used in 
the oxidation energy calculations of Eq. ([7]) . Compared to 
the tetravalent dioxide, A4O8, we observe that in A40g, 
the valency trend is from pentavalent in the early ac- 
tinides to tetravalent in the late actinides, indicating that 
the inclusion of O impurities into the dioxide leads to /- 
electron derealization in the early actinide oxides U40g, 
Np40g, and Pu 4 Og, whereas from An^Og onwards the 
compounds remain tetravalent. On the other hand the 
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A0,+ 1/8(V->A0 225 
AO -->A0, _+ 1/40, 



Am Cm 
Actinide ion 



FIG. 6: Oxidation (circles) and reduction (squares) energies 
of the actinide dioxides according to the definitions in Eqs. 
(J7J and ((9j). Energies are in eV per AO2 unit cell. 



fact that all the sesquioxides, A4O6, prefer the trivalent 
ground state configuration shows that the removal of O 
from the dioxide leads to localization of an additional 
/-electron. Whether the oxidation (Eq.©) or reduction 
(Eq.©) process actually takes place depends on the ther- 
modynamic conditions, in particular the O chemical po- 
tential. 

The results for the calculated oxidation and reduction 
energies are presented in Fig. [6l A binding energy of 
7.2 eV has been adopted for the O2 molecule (relative 
to the non spin-polarized O atoms. The atomic spin- 
polarization energy [gl] is 1.14 eV per O atom). The ac- 
curacy of the energy differences is not expected to be 
very high due to the ASA and the neglect of structural 
optimizations. Note, however, that for PU2O3, our cal- 
culated value of 2.6 eV/PuC>2 for the reduction process, 
Eq.©, lies within the range of LDA+U/GGA+U (U = 4 
eV) values calculated by Jomard et aZ. [llj] , 2.59/2.05 
cV/Pu0 2 , and Sun et al.^, 1.50/2.00 eV/Pu0 2 (data 
points taken from their figure). Nevertheless, here we will 
only concentrate on trends, which come out quite clearly. 

It emerges that E ox increases slightly and E red de- 
creases as we move through the actinide series, which 
reflects the increasing f-clcctron binding energy. In Ura- 
nium an f-electron is readily promoted into the valence 
band to facilitate uptake of super-stoichiometric oxygen, 
while this becomes increasingly more difficult with heav- 
ier actinides. On the other hand the U f-manifold is very 
reluctant to take on an extra f-electron to form the 3+ 
ion on reduction, while this becomes easier for the later 
actinides. The jump in reduction energy at Bk reflects 
the relative stability of the half-filled f 7 shell in BkOV 

According to the calculated energy differences, the only 
compound that oxidizes exothermically is UO2, whilst 
for all other actinide dioxides, the energy balance is pos- 
itive and increasing with actinide number. This result 
is in overall good agreement with experimental data. 
UO2+Z has been synthesized up to x=0.25 (L^Og). The 



crystal structure remains fluorite based, with the excess 
O atoms situated at interstitial sites. [631 With respect 
to Np, the only known stable binary oxide, apart from 
Np0 2 , is NP2O5, with a structure that can no longer 
be described as underlying fluorite. (70j In particular, the 
fluorite derived NP4O9 does not seem to exist. With 
respect to PuO^+z, the present calculations confirm the 
stability of PuC>2, in line with the majority of experimen- 
tal investigations on this material. The same conclusion 
was reached in an earlier SIC-LSD study [43| where the 
possible discovery of a higher composition binary oxide, 
Pu02+a; (x < 0.27) [7lJ was discussed. Beyond Pu there 
have been no reports of super-stoichiometric oxides. The 
positive slightly increasing oxidation energies in Fig. [5] 
confirm that the dioxide for all the trans-plutonium ac- 
tinides is increasingly stable against oxidation. 

Concerning the reduction reaction in Eq. (JHJ) , it is found 
that the energy balance always favors the dioxide. With 
respect to the early actinides U and Np the large value 
of the reduction energy is in agreement with the exper- 
imental fact that neither U2O3 nor NP2O3 exists. The 
sesquioxides from P112O3 to Cf2 03 do exist, however they 
are according to our calculations in principle unstable in 
air towards further uptake of oxygen and formation of 
the dioxide. In practice there may be appreciable barri- 
ers to the actual transformation, [72j which furthermore 
may be influenced by thermodynamic conditions (tem- 
perature and pressure) not considered here. 



IV. SUMMARIZING DISCUSSION 

From the study of the stable binary actinide oxides, 
a clear picture emerges that links the degree of oxida- 
tion to the degree of /-electron localization. In the early 
actinides, the /-electrons are less bound to the actinide 
ions which translates into valencies as high as 5+ and 
6+ for U-oxides for example. As one progresses through 
the actinide series, the /-electrons become increasingly 
tightly bound to the actinide ion, and eventually for Cf 
only the 3+ valency occurs naturally. The actinide ions 
play an active role in accomodating extra O, as their lo- 
calized /-electrons can act as electron reservoirs for the 
highly electronegative O-ion. In other words, whether 
the oxidation is favoured depends on the willingness of 
the actinide to delocalize further /-electrons. From the 
SIC-LSD study we find that the divalent configuration 
is never favoured, except maybe for EsO. Consequently, 
the monoxide never forms as the third electron is readily 
delocalizcd and made available for oxidation, which leads 
to the formation of trivalent sesquioxides. Whether the 
oxidation progresses further to produce the correspond- 
ing dioxide depends on derealization of the 4th electron 
on the actinide ion. As we have seen this derealiza- 
tion is less likely to happen for the late actinides, where 
the /-electrons tend to be more localized, and where ex- 
periment shows that the synthesis of Cf02 takes place 
under powerful oxidation conditions, e.g. in the presence 
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TABLE IV: Groundstate valency configurations of the ac- 
tinide and lanthanide oxides. Bold lettering indicates that 
the corresponding compound actually exists in nature. 
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of high pressure molecular or atomic oxygen. [72| For the 
early actinides, on the other hand, the oxidation to the 
dioxide occurs readily, and for example U2O3 and NP2O3 
do not exist naturally. Finally, the further oxidation from 
dioxide to higher oxide only occurs for UO2, where both 
the 5+ (U2O5) and 6+ (UO3) U-valencies exist, and for 
NpC>2 where the 5+ valency is obtained in NP2O5. We 
can speculate here as to whether the fact that we actu- 
ally determine a marginally pentavalent groundstate for 
UO2 (rather than the well established tctravalent config- 
uration) is an indication that an additional electron /- 
electron is on the brink of derealization, thus favouring 
the formation of higher U-oxidcs. The remaining dioxides 
are stable with respect to oxidation, as the gain in delo- 
calizing the 5th electron is no longer favorable compared 
to the corresponding loss in SIC energy. 

In Table IIVI we have collected our results of the SIC- 
LSD total energy calculations for the actinide oxides (the 
upper 4 rows). The numbers indicate the ground state 
valency configurations that we have determined for a 
given compound. Bold large letters are used to indi- 
cate those compounds that actually do exist in nature. 
It clearly emerges that only in those cases where the 
calculated ground state valency agrees with the nomi- 
nal charge expected from an ideal ionic picture, does the 
corresponding oxide seem to form. In other words, those 
oxides where our calculations predict a valency configura- 
tion that is not in agreement with simple charge count- 
ing will not form naturally, and the excess/shortage of 
charge indicates that oxidized/reduced compound will be 
favoured instead. These trends emphasize the very ionic 
nature of bonding in the actinide oxides. 

It is interesting to compare these trends in oxidation 
of the actinides to the corresponding behaviour of the 
lanthanides.[73l. |74T | In the latter the 4/ electrons arc 
overall more tightly bound to the lanthanide ion, which 
will favour lower oxidation numbers, compared to the 
spatially more extended 5/ electrons. The corresponding 
calculated valency configurations [73| are shown in Table 
IIVI (the lower 4 rows). The effect of the increased local- 



ization can be observed from the fact that a number of 
the lantanide monoxides actually exist, and it emerges 
from SIC-LSD studies that especially for the half-filled 
and filled shell /-electron systems, i.e. for EuO and 
YbO, the divalent configuration is energetically most fa- 
vorable. All the lanthanide sesquioxides occur in nature, 
and the ground state configuration of the corresponding 
lanthanide ions is trivalent. Even though the 4/-clcctrons 
tend to be more localized, the 3+ valency occurs natu- 
rally, which is in agreement with the trivalent configura- 
tion being the most favourable for the lanthanide met- 
als. With respect to further oxidation to the tetravalent 
dioxide, only Ce02, Pr02 (not shown in Table HV|). and 
TbC>2, are found to occur naturally or (for TbO^) under 
high O pressure, indicating that an additional /-electron 
only delocalizcs in the very early lanthanides, i.e. in Ce 
and Pr, where the /-electrons are less tightly bound, or 
in the middle of the series, where tetravalent TbC>2 is 
favoured by the half-filled shell configuration. [73j Higher 
oxidation numbers than IV do not exist for the lanthanide 
oxides. Again we observe that agreement between nomi- 
nal charge and calculated ground state valency is required 
for a given compound to exist. 



V. CONCLUSION 

We have studied the electronic structure of actinide 
oxides and specifically monoxides, sesquioxides and diox- 
ides, within the ab initio SIC-LSD band structure 
method. By studying the oxidation and reduction re- 
actions we have been able to conclude that the dioxides, 
from Np onwards, are the most stable compounds among 
the studied actinide oxides. Our study reveals a strong 
link between the preferred oxidation number and the de- 
gree of localization which is confirmed by comparing to 
the ground state configurations of the corresponding lan- 
thanide oxides. The ionic nature of the actinide oxides is 
reflected in that only those compounds can form where 
the calculated ground state valency agrees with the nom- 
inal valency expected from a simple charge counting. 
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